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ABSTRACT 

(Distribution  Limitation  Statement  A) 


The  effects  of  nuclear  scattering  and  ionization  energy  Toss  on  the  electro 
magnetic  fields  from  high-altitude  bursts  are  examined  u^ing  Monte  Carlo 
techniques  to  calculate  electron  trajectories.  Maxwell's  equations  and  the 
Lienard-Wiechert  potentials  are  both  used  to  evaluate  the  fields.  The 
obliquity  model  is  demonstrated  to  be  an  adequate  treatment  of  such  effects 
in  electromagnetic  pulse  calculations. 


AEWI.-TP-74-M 


CONTENTS 

Section 

* 

I  INTRODUCTION 

II  FIELD  EQUATIONS 

III  MONTE  CARLO  PROCEDURE 

IV  RESULTS 

REFERENCES 

DISTRIBUTION 


Page 

1 

2 

6 

9 

16 

17 


v 


! 


ArWL-7P-7«-45 


ILUJ'iTRATIOr.'S 

P igure  Page 

1  Electric  Field  Generated  at  the  Ground  by 

0.7-Mev  Electron  at  30  km  11 

2  Electric  Field  Generated  at  the  Ground  by 

0.53-Mev  Electron  at  30  km  12 

3  Electric  Field  Generated  at  the  Ground  by 

0.13-Mev  Electron  at  30  km  13 

4  [  lectric  Field  Generated  at  the  Ground  by 

0.7-Mev  Electron  at  20  km  M 

5  Electric  Field  Generated  at  the  Ground  by 

0.53-Mev  Electron  at  20  km  15 


AFWL-TR-74-46 


SECT  ion  I 
INTRODUCTION 

The  effects  of  nuclear  scattering  and  energy  loss  on  high-altitude  electro¬ 
magnetic  pulse  (FMP)  calculations  were  first  studied  by  G.  Knutson  (ref.  1)  in 
extensive  Monte  Carlo  scattering  calculations.  This  work  extends  the  Knutson 
work  by  calculating  in  a  simple  geometry  the  electric  radiation  fields  from 
Monte  Carlo  electron  calculations  similar  to  those  of  Knutson.  The  field 
calculations  are  done  simultaneously  using  the  high-frequency  approximation 
(ref.  2)  and  the  "particle”  formalism  (refs.  3,  4).  In  addition,  the  Mission 
Research  Corporation  (MRC)  obliquity  model  (ref.  5)  and  a  recent  modification 
are  examined.  Finally,  the  applicability  of  using  a  0.53-Mev  forward-scattered 
Compton  electron  to  simulate  the  effects  of  the  Compton  angular  distribution 
and  nuclear  scattering  for  incident  1.5-Mev  gammas  is  considered. 
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SECT  ion  II 
FIELD  EQUATIONS 


For  a  single  particle  of  charge  e  and  velocity,  v/c  =  B,  the  electric 
radiation  field  at  some  field  point  is  (ref.  6) 


E(r,t) 


e 

c 


,  r , 

n  x  L(n  -  B)  x  R 
R(1  -  n  •  B)' 


RET 


0) 


where  n  is  a  unit  vector  pointing  from  the  moving  charge  to  the  field  point, 

P  is  the  separation  between  the  particle  and  the  field  point,  and  [  Ippj 
means  that  the  enclosed  expression  is  to  be  evaluated  at  time  t  -  P/c  where  t 
is  the  time  at  the  field  point  of  interest.  For  an  initial  estimate  of 
scattering,  assume  a  horizontally  uniform  thin  layer  of  Compton  electrons 
deposited  at  30  km  at  t  =  0,  initially  moving  directly  downward.  If  there 
were  no  scattering,  the  electrons  would  rotate  about  the  magnetic  field  line 
with  angular  frequency  u  =  -eH/ym  c  where  H  is  the  ambient  magnetic  field  and 
Y  =  (1  -  v2/ c2)_l/2.  At  a  given  time  t,  then,  the  field  seen  at  the  ground 
is  found  by  summing  equation  (I)  over  all  electrons  in  a  disk  at  30  km  centered 
directly  over  the  field  point.  The  radius  of  the  disk  is 

r  =  [(3  x  10*  +  ct)?  -  (3  x  10fi)?]‘/?  cm 

Thus,  an  electron  at  the  center  of  the  disk  has  existed  the  longest,  i.e., 
t  seconds,  and  an  electron  at  the  edge  has  just  been  made  (Compton  scattered). 

Scattering  causes  the  electrons  to  deviate  from  the  circular  orbits  (ref.  7). 
This  degrades  the  coherence  cf  the  sources  and  results  in  a  lower  net  electric 
field. 

From  equation  (1),  the  transverse  radiation  field  from  a  particle  with 
velocity  B  =  v/c  is 
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ET(t>  =  e["x(("x  •  Gx)(By  •  eHxBz/nc) 

-  ("y-  By)  ex)  -  "z((V  By)  (6z  +  eWYmoc) 

■  (nz  ■  Gz)  (fty  "  eHx Vnc  ))  7(cR(1  -  n  • 

where  the  coordinate  system  has  the  z  axis  vertically  upward,  iix  is  the  earth's 

magnetic  field,  and  R  is  the  distance  from  the  source  particle  to  the  ground 

field  point.  ft  ,  and  ftz  are  the  accelerations  from  scattering  and  energy 

loss.  Due  to  the  specific  problem,  some  approximations  can  be  made.  R  is  taken 

as  30  km,  as  the  diameter  of  the  disk  containing  contributing  radiating  electrons 

is  always  small  (5  x  10"  cm  for  t  ■  t  -  R/c  <  2  x  10-8  sec)  compared  to  the 

altitude.  Also  1  -  n  •  B  »1  -  nzBz.  Finally,  nx  and  ny  can  be  taken  as  zero 

either  because  of  their  small  size  (n  ,  n  <  0.016)  or  because  the  integration 

a  y  — 

over  electrons  around  the  disk  results  in  zero  for  terms  in  equation  (2)  that. 

are  linear  in  n  and  n  .  Equation  (2)  can  now  be  rewritten, 
x  y 

Et(’p)  *  -*  [nz  (By  (»,  *  eHxey/YV) 

-  (°zBZ)  (6y  •eWYV))]/(cR('-nA)’)  (3) 

where  the  expression  on  the  right'is  evaluated  at  xp.  The  geometry  involved 
for  an  electron  at  some  point  on  the  disk  is  illustrated  by  the  following 
drawing,  a  is  the  altitude,  z  is  the  vertical  distance  traveled  by  the  electron, 
and  r  is  the  distance  from  the  center  of  the  disk  of  the  electron.  Assume  that 
the  electron  was  scattered  at  t  =  0.  The  proper  time  at  z  is  tp  =  t  +  z/c 
where  t  is  the  real  time  that  has  passed  in  the  life  of  the  electron.  The 
effect  of  the  particle  at  z  is  felt  at 

t  =  tp  +  (((a  *  z)!  +  -  (a  ♦  z)j/c  (4) 
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For  a  given  t  after  expanding  equation  (4) 

tp  ss  i  -  r?/2ac  (rj) 

D 

max 

rp(r)  Ej^ip^  dr  Ah  (6) 

where 

>*  ■  (?acT),/!  .  r.(r) 

is  the  density  of  Compton  electrons,  and  /.h  is  the  depth  of  the  deposition 
layer.  Ep(t)  is  the  net  field  on  the  ground  found  by  integrating  over  all 
contributing  electrons.  Now  rdr  -  -acdxp,  hence 

Ep(t)  =  27rac  /  pEj(xp)  dTp  Ah  (7) 

Jo 

The  field  calculation  is  based  on  the  high-frequency  approximation  (ref.  2). 
The  equation  for  the  transverse  field  is  for  a  given  proper  time  x. 


Ep(i)  =  2tt  f 
*'0 


4 


AFWL-TR- 74-46 


(8) 


where  J  is  the  current  and  is  expressed  in  terms  of  the  electron  motion  as 


J  '  ce  £  "lBy1  {(r  -  ri)/('  '  Bzl)  (9) 

i 

at  i,  is  the  particle  weight,  and  r^  is  the  particle  location  with  respect 
to  gamma  source.  Upon  integrating  equation  (8), 


(10) 


REp(f)  -  -2Tre  £  r.n.  By./(l  -  B2i)  (11) 

i 

EF(T)  .  -2ve  £  n,  By1(r,/R)/(l  -  B,,)  (12) 

i 

The  particle  calculation  is  for  an  infinite  source.  To  dc  the  field  calcula¬ 
tion  in  the  same  scenario,  the  gamma  source  is  taken  to  infinity.  Finally,  for 
fields 


Ef(t)  =  -27re  £  ni  By1/(l  -  Bzi) 
i 

and  for  particles 


Ep(x)  =  27iac  j  pEj(t)  dxp  Ah 


(13) 


(14) 


> 
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SECTION  III 

MOMTE  CARLO  PROCEDURE 

The  Monte  Carlo  calculation  of  the  motion  of  the  individual  particle 
includes  the  magnetic  field  turning,  elastic  nuclear  scattering,  and  the 
electron  energy  loss  from  ionization.  The  algorithm  is  almost  the  same  as 
used  by  Knutson  (ref.  1).  The  nuclear  cross  section  used  is 

do  =  2-nr r2  z2  (1  -  G2)  du/(c4(l  -  u  +  N)2)  (15) 

where  p  =  atom  density,  z  7.2,  u  =  cos  0,  rQ  =  2.82  x  10  1  cm,  and 

N  »  (1  -  B2)  (zl/3/(0.885)  ( 1 37 ) ) 2  (l.l3  +  3.76  (z/137B)2)/2P 

(16) 

The  energy  loss  formula  is 

dE/dx  =  A. 077  x  10"”  «e[ln(3.53  x  107  (y  -  1)(y?  -  1))  - 

(17) 

where  oe  is  the  ambient  electron  density. 

To  illustrate  the  Monte' Carlo  procedure,  look  at  a  particle  just  Compton- 
scattered  forward  with  velocity  Bz  and  y  s  (1  -  B2)  l^2.  The  total  nuclear 
cross  section  is  calculated. 

oT  -  1.04  x  10’22  p(l  -  B2)/(NB2)  (18) 

\ 

A  random  number  rn  is  calculated.  The  length  traveled  to  the  first  collision 
is 

d  =  -In (rn)/aT 
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The  real  time  passed  is 


At  ■  d/(3  x  10l°  B)  (20) 

The  z  coordinate  is  updated 

z  =  z  +  3  x  10,fJ  K  At/K  (21) 

where  K  is  the  obliquity  factor  developed  by  MRC  (ref.  5)  when  implemented. 

The  proper  time  is  also  updated 

7  =  t  +  z/(3  x  10‘°)  (22) 

Next  the  particle  is  turned  in  the  magnetic  field  for  a  time  t  and  the  contri¬ 
bution  to  the  integral  in  equation  (14)  is  calculated  for  the  particle.  Two 
more  random  numbers  n  and  r2  are  calculated  to  determine  the  scattering  angles 

cos  (o)  =  (2r,(l  +  N)  -  N ) / ( 2 r ,  +  N) 


0  =  2Tir2 

sin  <>  -  (1  -  cos7  d)'/'  (23) 

The  velocity  components  are  now  adjusted. 

a2  *  cos  ($)  sin  (6)/(b*  +  B*)1/2  (24) 

aj  «  B  sin  ($)  sin  (e)/(B*  +  B*),/2  (25) 

Bz  -  Bz  cos  (6)  -  a iBy  -  a2BxBz  (26) 

By  *  By  cos  (6)  +  a iBz  -  a2BxBy  (27) 

B>  =  Bx  cos  (e)  +  a2(B7  +  B7)  (28) 


) 
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Next,  the  energy  loss  is  calculated  and  the  total  velocity  changes  due  to 
scattering  and  energy  loss  are  used  to  update  the  particle  integral  in  equation 
(14).  Finally,  the  accumulated  particle  result  and  the  field  calculation 
result  are  accumulated  in  the  proper  time  bin  spanned  since  the  last  collision. 
The  MRf,  obliquity  factor  is  used  and  updated  in  the  field  calculation  when 
desired.  Each  subsequent  collision  of  the  particle  is  done  the  same  way  until 
2  x  l(Te  sec  of  proper  time  has  passed.  In  general ,  the  high-frequency 
approximation  and  the  "particle"  model  agreed  to  5  percent. 
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SECTION  IV 
RESULTS 

Figures  1  through  3  demonstrate  the  results  at  30  km  for  various  initial 
energies.  Figure  1  is  for  0.7  Mev  and  is  the  most  similar  to  the  typical 
energy  used  in  more  sophisticated  codes.  The  circular  orbit  result  at  2  x  10"8 
sec  is  about  2.9  times  larger  than  the  Monte  Carlo  scatter,  and  the  obliquity 
results.  The  agreement  between  the  Monte  Carlo  and  the  obliquity  is  excellent. 
Figure  2  is  for  0.53  Mev,  and  the  decrease  due  to  scattering  from  circular 
orbits  is  a  factor  of  0.26.  The  agreement  between  the  scatter  and  the  obliquity 
is  not  as  good  as  for  higher  energies  but  is  still  quite  acceptable.  The  square 
on  figure  2  is  a  result  by  Sollfrey  (ref.  8)  for  the  same  case.  Finally,  figure 
3  is  for  0.13  Mev.  The  fractional  decrease  from  circular  orbits  is  about  0.2. 
The  discrepancy  between  scatter  and  obliquity  is  now  more  pronounced. 

Figures  3  and  4  contain  the  same  results  as  figures  1  and  2  at  20  km.  The 
agreement  between  the  obliquity  factor  and  the  Monte  Carlo  is  still  reasonably 
gocd. 

In  addition  to  comparisons  of  obliquity  to  Monte  Carlo  scatter,  other 
in/estigations  were  made.  At  a  1  November  1973  meeting  at  AFWL,  Conrad 
Longmire  of  MRC  described  a  new  obliquity  model.  If  the  equation  of  change  of 
the  old  obliquity  factor  is  represented  as 

dN/dt  =  F ( E , o ) 

where  E  is  the  particle  energy  and  p  is  the  density,  then  the  new  equation  is 

dN/dt  =  NF(E,p) 

Results  are  plotted  in  figures  2,  4,  and  5.  Results  were  not  plotted  on 
figure  1,  as  the  new  model  did  not  change  significantly  from  the  old.  In  all 
cases,  the  new  model  is  a  significant  improvement  over  the  old.  The  effect 
of  using  the  new  model  as  compared  to  the  old  was  tested  in  a  full  EMP  calcula¬ 
tion  with  CHEMP,  the  AFWL/DYT  EMP  code,  and  only  very  slight  differences  were 
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noted  in  the  tails  of  the  EM  pulses.  This  means  that  older  calculations  with 
the  old  model  are  still  useful  and  also  that  field  results  are  not  very 
sensitive  to  fine  details  in  the  source. 

Plotted  in  figures  2  and  5  are  the  results  of  combining  the  Klein-flishina 
angular  distribution  from  an  incident  1.5-Mev  garrma  with  the  Monte  Carlo 
scattering  and  with  the  new  obliquity  model.  This  was  done  to  test  the 
equivalence,  as  stated  by  Cullen  Crain  {ref.  4),  between  the  Klein-Nishina 
plus  scatter  and  a  0.53-Mev  forward-directed  electron  with  scatter  alone.  At 
3^  km  the  agreement  is  very  good.  At  20  km  at  2  x  10’-  sec,  however,  the 
0.53-Mev  electron  is  about  30  percent  below  the  comparison  result.  Thus,  the 
use  of  the  0.53-Mev  forward-directed  electron  to  model  the  currents  from  a 
1.5-Mev  gamma  pulse  would  appear  to  underestimate  the  fields. 

In  summary,  the  obliqqity  models  do  an  adequate  job  in  simulating  the 
effects  of  nuclear  scattering  with  the  new  model  being  clearly  better  than  the 
old.  These  results  do  not  establish  the  obliquity  result  over  all  possible 
regions,  but  they  do  substantiate  the  much  more  extensive  Monte  Carlo  results 
by  Knutson  (ref.  1). 
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Figure  1.  Electric  Field  Generated  at  the  Ground  by  0.7-Mev  Electron  at  30  km 
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Figure  2.  Electric  Field  Generated  at  the  Ground  by  0.53-Mev  Electron  at  30  km 
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Electric  Field  Generated  at  the  Ground  by  0.13-Mev  Electron  at  30  km 
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Figure  4.  Electric  Field  Generated  at  the  Ground  by  0.7-Mev  Electron  at  20  km 
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Figure  5.  Electric  Field  Generated  at  the  Ground  by  0.53-Mev  Electron  at  20  km 
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The  effects  of  nuclear  scattering  and  ionization  energy  loss  on  the  electromagnetic 
fields  from  high-altitude  bursts  are  examined  using  Monte  Carlo  techniques  to 
calculate  electron  trajectories.  Maxwell's  equations  and  the  Llenard-Wlechert 
potentials  are  both  used  to  evaluate  the  fields.  The  obliquity  model  Is  demon¬ 
strated  to  be  an  adequate  treatment  of  such  effects  In  electromagnetic  pulse 
calculations. 
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